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Abstract 

The phase structure of four-fermion theories is thoroughly investigated with 
varying temperature and chemical potential for arbitrary space-time dimensions 
(2 < D < 4) by using the 1/N expansion method. It is shown that the chi- 
ral symmetry is restored in the theory under consideration for sufficiently high 
temperature and/or chemical potential. The critical line dividing the symmetric 
and broken phase is given explicitly. It is found that for space-time dimension 
2 < D < 3 both the first-order and second-order phase transition occur depending 
on the value of temperature and chemical potential while for 3 < D < 4 only the 
second-order phase transition exists. 
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1 INTRODUCTION 



The idea of spontaneous symmetry breaking has played a decisive role in recent devel- 
opments of particle physics. In fact the spontaneous breaking of the gauge symmetry 
SU (3) c <8> SU (2)2, <8> U (l)y is the basic ingredient of the standard electroweak theory and 
grand unified theories (GUT) are constructed on the basis of the Higgs mechanism. In 
these approaches, however, the spontaneous breaking of symmetries is introduced phe- 
nomenologically in terms of the Higgs field. In other words the search for the dynamical 
origin of the spontaneous symmetry breaking is passed over by introducing elementary 
Higgs fields. Under this circumstance it is very interesting to look for a possible dynamical 
mechanism of the spontaneous symmetry breaking. 

More than 30 years ago Y. Nambu and G. Jona-Lasinio first introduced the idea 
of dynamical symmetry breaking in particle physics^ and the idea has attracted a vast 
number of researchers in particle physics. One of the major applications of this idea has 
been made on the construction of composite Higgs models such as the technicolor model.! 
It is important to test phenomenologically whether any of these composite Higgs models 
is a candidate of the underlying theory of the standard electroweak theory. 

One of the possible environments where the composite Higgs models may be tested is 
found in the early universe where the symmetry of the primary unified theory is broken 
down to yield lower-level theories which describe phenomena at lower energy scales. In the 
early universe it is not adequate to neglect the effect of the curvature, temperature and 
density. Hence it is important to examine a theory in the finite curvature, temperature 
and density. In our study we are interested in the spontaneous breaking of the electroweak 
symmetry and hence we consider the scale much smaller than the Planck scale. In such a 
situation the essential feature of the early universe may be described only with the effects 
of finite temperature and density, the curvature effect being negligible. In the present 
communication we would like to report our investigation in a simple model field theory 
for composite Higgs fields at finite temperature and chemical potential. 

We work in a four-fermion theory in the flat space-time with arbitrary dimensions. 
The theory is composed of iV-component fermions with four-fermion interactions. We 
employ the 1/N expansion method to estimate the effective potential for the fermion- 
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antifermion composite field. In two space-time dimensions the theory reduces to the 
Gross-Neveu modeli which possesses the discrete chiral symmetry that is broken dy- 
namically. In three space-time dimensions the theory defines a model which is solvable 
and renormalizable in the sense of the 1/N expansion.^ Beyond four space-time dimen- 
sions the theory is nonrenormalizable. We confine ourselves to the space-time dimensions 
greater or equal to 2 and less than 4. 

We introduce the temperature and chemical potential in the theory through the stan- 
dard procedure. We investigate the behavior of the effective potential in the leading order 
of the 1/N expansion by varying the temperature and chemical potential for arbitrary 
space-time dimensions. Through the study of the shape of the effective potential and the 
detailed analysis of the gap equation we observe the phase transition from the phase with 
the broken chiral symmetry to the phase with the chiral symmetry when the temperature 
and chemical potential vary. With this analysis we shall be able to derive explicitly the 
critical curve on the temperature-chemical potential plane. We find that for space-time 
dimensions less than 3 both the first-order and second-order phase transition exist while 
for space-time dimensions greater or equal to 3 only the second-order phase transition is 
realized. At some specific points on the critical curve we obtain an explicit expression 
for the critical temperature and/or chemical potential. We also calculate the dynamical 
fermion mass as a function of the temperature and chemical potential. 

In section 2 we briefly review the general properties of four-fermion theories in space- 
time dimensions 2 < D < 4 with vanishing temperature and chemical potential. In 
section 3 we introduce the temperature and chemical potential in the theory. We calculate 
the effective potential in the leading order of the 1/N expansion and observe the change 
of the phase. We derive the critical curve which divide the chiral symmetric phase and 
the chiral asymmetric phase. The section 4 is devoted for the concluding remarks. 

2 FOUR-FERMION THEORY IN ARBITRARY SPACE-TIME 
DIMENSIONS 

In the present section we briefly summarize the characteristic features of four-fermion 
theories in arbitrary space-time dimensions 2 < D < 4 for vanishing temperature and 
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chemical potential. We start with the simple Lagrangian 

N \ N 

£ = E favrvWk - -£ £(^) 2 , (2.1) 
k=i ZiV k=i 

where index k represents the flavors of the fermion field ip, N is the number of flavors 
and Ao is a bare coupling constant. In the following, for simplicity, we neglect the flavor 
index. In two dimensions the theory is nothing but the Gross-Neveu model! in which 
Lagrangian Q2.1] ) is invariant under the discrete chiral transformation, 

^-^75^- (2.2) 
The above discrete chiral symmetry in two dimensions prevents the Lagrangian to have 



mass terms. In arbitrary dimensions the transformation (2.2) may be generalized so that 

ipijj — > —i/ji/j. 

The theory has a global SU(N) flavor symmetry : the Lagrangian is invariant under 
the transformation, 

,j, — > e iS ^ T > (2.3) 

where T a are generators of the SU(N) symmetry. Under the circumstance of the global 
SU(N) symmetry we may work in the scheme of the 1/N expansion. 

For practical calculations in four-fermion theories it is convenient to introduce auxil- 
iary field a and consider the following equivalent Lagrangian 

C y = jk-y^ip - ^-o 2 - ^ . (2.4) 
zAo 

If the non-vanishing vacuum expectation value is assigned to the auxiliary field a, then 
there appears a mass term for the fermion field ip and the discrete chiral symmetry (the 
Z2 symmetry in odd dimensions) is eventually broken. 

We would like to find a ground state of the system described by four-fermion theories. 
For this purpose we evaluate an effective potential for composite field ifjifj in the theory 
described by Eq. ( |2.1| ). This effective potential is essentially the same as the one for field 
a in the theory described by Eq. ( |2.4| ). The effective potential Vo(cr) in the leading order 
of the 1/N expansion reads (See Appendix A) 

V (a) = ^-a 2 + zln det(^y9 M - <r) + 0(1/JV) , (2.5) 



where the suffix for Vo(cx) is introduced to keep the memory that T = fi = with T 
the temperature and \x the chemical potential. 

For integral dimensions the effective potential is in general divergent. Performing the 
integral in the potential (|2.5[) we obtain 



V M = _L a 2 _ 1 r ( 1 - - )a D . (2.6) 
V ' 2A {2k) d / 2 D \ 2 J K ' 

We clearly see that the effective potential is divergent in two and four dimensions. It 

happens to be finite in three dimensions in the leading order of the 1/N expansion. If 

the next-to-leading order is taken into account, it may be divergent in three dimensions 

D = 3. 

As is well-known, four-fermion theory is renormalizable in two dimensions. Therefore 
the potential (|2.6|) is made finite at D = 2 by the usual renormalization procedure. 
For D = 3 four-fermion theory is known to be renormalizable in the sense of the 1/N 
expansion.! The potential (|2.6| ) is finite by itself and hence we do not need renormalization 
in the leading order of the 1/N expansion at D = 3. In four dimensions four-fermion 
theory is not renormalizable and the finite effective potential can not be defined in four 
dimensions. We regard the effective potential for D = 4 — e with e sufficiently small 
positive as a regularization of the one in four dimensions. 

We perform renormalization in two dimensions by imposing the renormalization con- 
dition, 



d 2 V n (a) 



da 2 

where o"o is the renormalization scale. From this renormalization condition we get 



' ' ' p-i)r(i-^K~ 2 - (2.* 



A A r (2tt) d / 2V V 2 

The renormalized effective potential reads 
VM 1 a 2 



2Xa 2 



h^P-i)rfi-?K-7^rfi-?K. (2.9) 



2(2tt)^/ 2V ' V 2) al (2n) D / 2 D \ 2 J o$ 
Here we used the dimensionless coupling constant A defined by 

\ = \ r a°- 2 . (2.10) 
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The renormalized potential ( |2.9| ) is no longer divergent in the whole range of D considered 
here : 2 < D < 4. 

The ground state of the theory is determined by observing the minimum of the effec- 
tive potential. The necessary condition for the minimum is given by 



dV (a) 



da 



ma, 



,D-2 




1 



I ( 1 — — I 

2 / /m N ■ - 



A A c (2tt)^/2 Wo 



0. 



(2-11) 



where m is the dynamical fermion mass and A c is defined by 

A c =(2vr)^ 2 \(1-D)T (l-y ^ 



(2.12) 



If the coupling constant A is no less than a critical value A c , the gap equation allows a 
non-trivial solution which is given by 



m = ctq 



(27T) g / 2 



l/(2-D) 



(2.13) 



In Fig. |l|. we plot the critical coupling constant A c as a function of dimension D . The 
non-trivial solution m of the gap equation fl2.11|) corresponds to the nonvanishing vacuum 
expectation value of the composite field a and is equal to the dynamically generated 
fermion mass. For some special values of D the solution m simplifies : 



m = a e 



m — Co ( 2 



1-tt/A 



D = 2. 



:D = 3 



(2.14) 
(2.15) 



and A c also simplifies : 



A c = 
A - n 

c "71 



;D = 2, 
:D = 3. 



(2.16) 
(2.17) 



As is well-known, the shape of the effective potential Vo(er) is of single-well for A < A c 
while it is of double-well for A > A r . 
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Fig. 1 Critical coupling A c of the four-fermion theory as a function of dimension D . 



The renormalization group f3 function is found to 

dX 
(T — 
aao 

D-2 



X, 



A(A C - A) 



(2.18) 
(2.19) 



in the leading order of the 1/N expansion. It simplifies for some special values of D : 

7T 

;D = 3. 



/9(A) 



-A - l 



7T 



(2.20) 
(2.21) 



We find from Eq. ( ^.19| ) that A = A c is the ultraviolet stable fixed point. In particular 
for D = 2 the theory is asymptotically free and for D = 3 the theory has a nontrivial 

7T 

ultraviolet fixed point at A = —j= . 

v2 

We turn our attention to the theory governed by the following Lagrangian, 

C = - T^K^) 2 + (^75^) 2 ] ■ (2.22) 

In four dimensions this Lagrangian defines the Nambu-Jona-Lasinio model! The La- 
grangian is invariant under the chiral U(l) transformation in even dimensions, 



(2.23) 



The chiral U(l) symmetry prevents the Lagrangian to have mass terms. Using the 
auxiliary field method, we obtain the effective potential 

V M = ^ a ' 2 + tln det (*7^ - °') + 0(l/N) , (2.24) 
2Aq 



a 



V(J 2 + vr 2 . (2.25) 



Corresponding to the two four-fermion interaction terms in the Lagrangian ( j2.22j ), we 



need two kinds of auxiliary fields a and ir. In two dimensions it is well-known that it 
loop has an infra-red divergence and we can not neglect the next to leading order terms 
of the 1/N expansion.! Therefore our present analysis is not sufficient to deal with the 
Lagrangian fl2.22j ) in two dimensions. 

The effective potentials (|2.5| ) and ( |2.24j) have the same form in the leading order of 



the 1/N expansion. Thus using the effective potential (|2.5|) , we can discuss the phase 
structure of both theories except for the case in two dimensions. 

3 PHASE STRUCTURE AT FINITE TEMPERATURE AND 
CHEMICAL POTENTIAL 

3.1 Effective potential and gap equation 

As we have seen in the previous section the chiral symmetry is broken above the critical 
coupling. The symmetry broken spontaneously may be restored in an environment of 
high temperature and density. To see whether this situation is realized we would like to 
study the phase structure at finite temperature and chemical potential. 

The n-point Green function at finite temperature and chemical potential is defined 

by 

r0 , _ E^e-^-^VroOrQ ■■■ <f>(x n ))\a) 

n ^ e -P(E a -(iN a ) ' 

where <f>(x) represents a component of fermion field ip{x) or ip(x), E a and N a are the 
energy and particle number in the state specified by quantum number a respectively, 
(3 = - — with k the Boltzmann Constant and T the temperature and \x is the chemical 

rC J. 



potential. The generating functional for Green functions (|3.1| ) reads 

(3.2) 



Z^[J] = Tr e-^-^Texpi J dxJ(x)(f)(x) 
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where H is the Hamiltonian, iV is the number operator and J(x) is the source function 
for field (j>{x). 

Following the standard prooednre of the Matsnbara Green functkj (See Appendix 
B) we calculate the effective potential in our theory in the leading order of the 1/N 
expansion as described already in the case of T = fi = in Eq. fl2.5| ). We find 

1 9 lr-, / d D ~ 1 k , H — a , 

Ho) = n/ + 1 pVfw)^ lDdet —- (3 ' 3) 

where the four-momentum k^ is given by 

hP = (iu n -n,k), (3.4) 
and the discrete variable uo n given by 

Un = ~p , (3.5) 

according to the anti-periodic boundary condition. 

Performing a summation and integrating over angle variables (See Appendix B) we 

get 

V(a) = V (a) + Vto(c) , (3.6) 

where Vo(cr) is the effective potential for T = \i = shown in Eq. ( |2.6|) and V^^{a) is 
given by 



V^{a) 



1 y/2 1 D 



^( 27r )(D-D/2 r ^-l x ff o 



x / dkk D ~ 2 { In— ^r-^- + In — m ; } . (3.7) 

It is found that a) is finite in the space-time dimension 2 < D < A. Thus we need not 
renormalize the parts including the effects of finite temperature and density. We apply 
the same renormalization condition as shown in Eq. (2.7) and obtain the renormalized 



effective potential as in the T = \x = case. The renormalized effective potential V(a) 
is given by replacing V in Eq. ( |3.6|) by the one in Eq. ( |2.9|) . 
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If we perform the integration first in Eq. ( |3.3| ) and leave the infinite sum we find 

1-D" 



2-D/2-1 r 



2A (3 (4tt)( d - 1 )/ 2 

oo 

E [{K + + a 2 } (D - 1)/2 - {K + ^) 2 } (z? - 1)/2 



(3.f 



Although the above expression is essentially equivalent to the one in Eq. (|3.7|) , we shall 
often use the expression ( p.8|) in the later argument of the gap equation. 
The gap equation is obtained through the stationary condition 

dV{a) 



da 



0. 



(3.9) 



where is a possible dynamical mass of the fermion. If we apply the expression 
in Eq. (|3.9| ), we find 

A A c (27t)^/ 2 I 2/U 
s/2 



2-D 



T((D- l)/2)(27r)^- 1 )/ 2(Jo 
1 f 



x / dfc/c 

JO 



D-2 



+ 



0. 



(3.10) 



In Eq. ( |3.10| ) we neglected the trivial solution = 0. If Eq. ( |3.8|) is employed, the 
following expression is obtained, 

{3 — D\ 



1 1 

a "a: 



(3a (47t)^-D/ 2 n 



E 



4 



0. (3.11) 



3.2 Phase boundary 



The expression ( |3.7| ) is useful to calculate the effective potential numerically as a function 
of the composite field a . Using Eq. ( |3.7| ) we calculate the effective potential for D = 2 
to find the behavior of the effective potential! as shown in Fig. 0. 
As is seen in Fig. |2|, we observe that 
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0.2 0.4 0.6 0. 
kT/m 

(a) The T-/J, plane showing 
the fixed /i direction A 
and the fixed T direction B. 



V/m 2 



0.004 - 



0.002 



-0.002 



-0.004 - 








0.2 



0.4 0.6 

a/m 

(c) Behavior of the effective 
potential as varies 
in the fixed T direction B. 



0.8 



V/m 2 



0.004 - 



0.002 



-0.002 



-0.004 - 




0.2 0.4 0.6 0.8 

a/m 

(b) Behavior of the effective 
potential as T varies 
in the fixed fi direction A. 




0.4 0.6 

fi/m 

(d) The dynamical fermion mass 
as a function 
of /i with T fixed. 



Fig. 2 Behaviors of the effective potential and the dynamical fermion mass for D = 2 
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1. There is a second-order phase transition as T is increased with jj kept small fixed. 

2. There is a first-order phase transition as [i is increased with T kept small fixed. 

Thus we expect that the first- and second-order phase transition coexist on the T — // 
plane. For D just above 2 this property is expected to persist. On the other hand at 
D = 3 we observe only the second-order phase transition for varying T and /!.! This 
situation is illustrated in Fig. ||. 



V/m 3 




-0.002 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

a/m 

Fig. 3 Behavior of the effective potential at D = 3 along the B line in Fig. 2 (a). 

The above preliminary analysis suggests that the first-order phase transition observed 
at D = 2 disappears at some critical dimension D c < 3. 

In order to see the situation more precisely we would like to perform a rigorous 
analysis on the critical temperature T c and chemical potential fi c . We would like to find 
the critical line on the T — // plane in an analytic form if possible. Let us first assume that 
the phase transition is of the second order. Since for the second-order phase transition 
the nontrivial dynamical fermion mass develops continuously at the critical temperature 
and chemical potential as illustrated in Fig. 4, the critical temperature and chemical 
potential are determined by letting — > in the gap equation. Thus we set = 
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Fig. 4 Behavior of the dynamical fermion mass in the case of the second-order 
phase transition as T varies with \x fixed . 



in Eq. (|3.11| ) and obtain 



1 1 

a ~ r 



(3.12) 



where we used the following relation 
with the definition of the generalized zeta function ((z, a),0 



S \ ' 2 2vr / S I ' 2 2vr 



(3.13) 



C(«,a) = 



n=0 



n + a) 



Re z > 1 



(3.14) 



Eq. (|3~lp provid es us with the relation between the critical temperature T c and chemical 
potential fi c , i. e. the critical curve on the T — fi plane if the transition is of the second- 
order. It should be noted that the relation ( p,13|) holds only for D < 2 according to the 
condition Re (2v) > 1 as given in Eq. (|3.14 ). In Eq. ( |3.12|) , however, we analytically 
continue the variable D to the region KeD > 2. As a typical example the second-order 
critical curve at D = 2.5 is presented in Fig. [5] using Eq. ( |3.12| ). 
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fi/m 
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n 1 1 r 



j i i i i u L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0. 
kT/m 

(a) Second-order critical curve 
at D = 2.5 
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d=3.2 
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_D=2^ 




i i 


1 i II 1 li 







0.2 



0.8 



0.4 0.6 

kT/m 

(b) Second-order critical curves 
at D = 2.0,2.5,2.8,3.0,3.2,3.5 



Fig. 5 Second-order critical curves 



Eq. (|3.12| ) simplifies when D = 2 and D = 3. We note that the generalized zeta 
function a) has a pole at z — 1 and its behavior near z — 1 is given by 

f 



CO, «) 



ip(a) 



(3.15) 



where ip(a) is the digamma function. Using Eq. ( |3.15| ) for z = 3 — D and applying some 
formulae for the digamma function we find for D = 2 



2Re^(l + ^)-Re^(l + ^)=ln^ 



(3.16) 



where we employed Eq. ( p. 13 ) with m the dynamical fermion mass for T = /i = 0. Eq. 
( |3.16|) exactly reproduces the result obtained in Ref. 11. To derive the corresponding 
formula for D = 3 we employ the following expansion formula for the generalized zeta 
function: 

((x, a) = - - a + x jln T{a) - ^ In 2vr| + 0{x 2 ) . (3.17) 
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Using Eq. ( [3.171 ) with x = 3 — D and a 



+ i-J^- we rewrite Eq. (|3.12| ) as follows, 



2tt 



1 

A 



1 



V2 



In 



./V 



TV 



r 



./V 

' j 2tt 



(3.18) 



which reduces to the simple formula 



a /3m/2 



2 cosh — . 
2 



(3.19) 



Eq. ( |3.19| ) exactly agrees with the one obtained in Ref. 9. 

Our formula ( |3.12| ) gives a correct critical curve on the T — fi plane as far as the 
transition is of the second order. If the transition is of the first order, the curve given by 
the formula fl3.12|) corresponds to the local maximum (first extremum) of the effective 
potential and does not give a phase boundary. In such a case we have to find the 
critical curve corresponding to the true minimum (second extremum) of the effective 
potential. The critical curve for the first-order phase transition is obtained by eliminating 
the variable a in the following two equations: 

dV{a) 



V(a) = 



da 



0. 



(3.20) 



where we assume that a^O. 

As we have seen before there exists a first-order phase transition for 2 < D < 3. Here 

we have to perform an analysis based on Eq. ( |3.20| ). Unfortunately we are unable to 

solve this problem analytically. We then employ expression ( |3.6| ) with ( |3.7| ) and ( |3.10| ) 

to find the critical T and /i numerically through the use of MATHEMATICA on our 

workstation. For example by numerically calculating the effective potential at D = 2.5 

for (3m = 10 we obtain a figure of the effective potential as shown in Fig. 6. 

As can be seen in Fig. 6 the critical chemical potential is given by — = 0.81 at 

m 



— = 0.96. Repeating this kind of analysis many times we find the critical curve for 
m 

the first-order phase transition at D — 2.5. In Fig. 7 is shown the first-order critical 
curve by the dashed line while the second-order critical curve is given in the full line. We 
performed the full analysis of the above type for 2 < D < 3 and obtained the critical 
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V/m 



2.5 



0.0025 - 



0.002 - 



0.0015 - 



0.001 - 



0.0005 




0.2 0.4 0.6 

a/m 



0.8 



Fig. 6 Searching for the first order critical point by directly observing the behavior of 
the effective potential. 



fi/m 




0.4 
kT/m 

Fig. 7 Critical curve in D = 2.5. The dashed line represents the first-order phase 
transition while the full line represents the second-order phase transition line. The dotted 
line signals an appearance of the first extremum (local maximum) in the effective potential 
and has nothing to do with the phase transition. 
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9 e — — — _ ^\ 












i i 


Iill 1 1 1 



0.2 0.4 0.6 0.8 1 

kT/m 

Fig. 8 Critical curves for 2 < D < 4 . 

curves as shown in Fig. |[ In Fig. |8| we clearly observe that the first-order critical line 
given in the dashed line smoothly disappears at D = 3. We in fact checked that the 
phase transition is always of second order for 3 < D < 4. 



3.3 Specific points on the critical curve 

It is interesting to note that we are able to find analytically some specific points on the 
critical curve. Those points are shown in Fig. 7 by A, B, C and D. 

Let us first consider point A. Since we know that at this point the transition is of 
second order, we are free to use the formula ( P-12| ) for determining point A. We set \x = 
in Eq. ( |3.12j ) to obtain 



1 1 
A ~ AT 



2(2^/- r ^ ( h tl f "1^ ' 



1)C(3-D) = 



(3.21) 



fir V 2 

where ((z) is the zeta function. Using Eq. ( p,13|) we rewrite Eq. ( p.21|) in the following 
form: 

2r 



pm = 2tt 



3-D 



D 



(2^-l)C(3-D) 



l/(D-2) 



(3.22) 



Eq. ( |3.22D gives the critical temperature. 
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For D = 2 Eq. ( |3.22| ) reduces to the well-known formula for the critical temperaturejll 



(3m = ne 



(3.23) 



To see this we employ the following expansion formulae in Eq. (|3.22|) 

r(e) = i- 7 + 0(e 2 ), 



rQ + e) = v^{l-£(21n2 + 7 ) + 0(£ 2 )} 



C(l + 2e) = -+7 + 0(e), 
with 7 the Euler constant and to obtain 

(3m = lim2vr {1 + 2e(ln2 + 7)}~ 1/2e 

£— >0 



(3.24) 



(3.25) 



which clearly reproduces Eq. (|3.23|) . For D = 3 Eq. (|3.22|) again reproduces the well- 
known formula! 



(3m = 21n2. 



(3.26) 



To derive Eq. ( |3.26| ) from Eq. Q3.22| ) we only need the following expansion formula near 
D = 3, 



C(3 -D) = - 1 -- 1(3 - £>) ln27T + 0((3 - D) 2 ) 



(3.27) 



In Fig. |^ the critical temperature is plotted as a function of dimension D by the use of 
Eq. (fE22p. 

We next consider point B. On the /i-axis for which T = {(3 — > oo) the gap equation 
(|3.10|) reads for jjl < 

2- 



(3.28) 



and for // > m^ 
T .'2-D 



D-2 l>-2, -\'7 /V"" 



D - 1 



UD-2 

dk ; = = . (3.29) 



fc 2 + 



Eq. ( |3.28|) gives a simple solution 



(3.30) 
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Fig. 9 critical temperature for ji = as a function of dimension D 
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Fig. 10 Dynamical fermion mass as a function of /x with T = 
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which is typically of the first-order type since it has a gap at /j = fi c as is seen in Fig. 
0(d). On the other hand Eq. ( |3.29| ) gives a more complicated solution which is shown in 
Fig. [K| This solution is typically of the second-order type. 

For 2 < D < 3 the phase transition along the /i-axis is of first order as we have 
seen before. Hence the solution ( |3.30| ) gives a true vacuum and the solution of Eq. ( |3.29| ) 



corresponds to the first extremum of the effective potential. We adopt the solution (|3.30|) 
to study the condition, 

lim V(a = m) = 0. (3.31) 



Using Eq. ( ^.6[ ) with Eq. fl3.7p we rewrite the condition ( |3.31| ) as follows, 



1 r ( i— -) m D + r dkk°- 2 (k -fi) = 0, (3.32) 



_l_ 1 I lib , 7— » 

D V 2 J T ( D 



o 



2 

which provides us with the critical chemical potential for 2 < D < 3 , 

ll = m lh(^R,£±l)\ 1,D . (3,3) 



.4 V 2 2 

It is easy to check that Eq. ( |3.33| ) satisfies the condition . For D = 2 Eq. ( |3.33| ) simplifies 
toiJH 

m 

M — ~^j= ■ (3-34) 
For 3 < D < 4 the phase transition along the //-axis is of second order and hence we 
have to adopt the case of Eq. ( |3.29| ). Letting — > in Eq. ( |3.29p and solving for /x 
we find 

B . (3.35) 



.2 V 2 2 

Obviously Eq. ( |3.35| ) satisfies the condition \i > mp^. For D = 3 Eq. ( |3.35| ) simplifies to 



H = m , (3.36) 

which agrees with the known result .i In Fig. [TT] the the critical chemical potential is 
plotted as a function of dimension D by using Eqs. (|3.33|) and ( |3.35| ) . 



It is worth noting here that, if we would have used the second- order result ( p.35|) for 



2 < D < 3 , we would have obtained point D shown in Fig. 7. The value of \i at point 
D is given by Eq. (|3~35p with 2 < D < 3 . For D = 2 we find 

A* = y > (3-37) 
20 



fi c /m 




Fig. 11 Critical chemical potential for T = as a function of dimension D . 

instead of Eq. (|3.34| ). This, however, corresponds to the first extremum of the effective 
potential and has nothing to do with the phase transition.0 

Finally we consider point C on the critical curve. This point appears only when 
2 < D < 3 . At point C the first-order critical curve meets with the second-order critical 
curve. To find out the value of T and fi at point C we proceed as follows: Let us consider 
the gap equation 



dV_ 



0. 



(3.38) 



We fix the dynamical mass at the first nonvanishing extremum . Thus we have 



m=o- 



(3.39) 



By fixing a at o m temperature T is related to chemical potential fi through Eq. (|3.39|) . 
As shown in Fig. 12 the relation between T and \i is plotted in the form of the curve 
G just as we have done for the critical curve. Since the dynamical fermion mass is fixed 
to cr m on this curve G, the behavior of the effective potential along the curve G looks as 
in Fig. 13. Obviously the derivative of the effective potential, i. e. f(a) , changes sign 
when the curve G crosses the critical curve. On the critical curve the effective potential 
is flat in the region a < a m as seen in Fig. 13. Thus at the crossing E of the curve G 
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Fig. 12 The T-/x plot as given by 
Eq. (3.39) when a is fixed at a m . 
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Fig. 13 Behavior of the effective 
potential as T and fi vary on the 
curve G [a m = 0.3). 



and the critical curve the second derivative of the effective potential has to vanish, 



d " V f( ^ °f n 
^ = /((T) + 0r ^ = °- 



Combining Eq. ( 3.39j) with Eq. ( 3.40| ) we have a condition for the point E, 

df 



da 



0. 



(3.40) 



(3.41) 



The point E divides the second-order region from the first-order region along the curve 
G. By letting a — > in Eq. ( |3.40| ) we reach to the point C in Fig. 7. Hence the condition 
to get the point C is given by 

0. (3.42) 



da 



<T=0 



Inserting Eq. (|3.11|) into Eq. (|3.42|) we finally obtain the condition for point C, 

0. (3.43) 
By solving Eq. (|3.43|) for we have a curve in the T-/i plane. Looking for the crossing 



of this curve with the critical line we find the value of (3 and /i at point C. In Fig. |1J the 
value of (3fj, at point C is plotted as a function of dimension D by using Eq. ( |3.43j ). 
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Fig. 14 The value of kT / \i at point C as a function of dimension D . 
3.4 Dynamical fermion mass at finite T and /i 

We have seen that the broken chiral symmetry at low temperature and chemical poten- 
tial is restored at certain critical temperature and chemical potential. The symmetry 
restoration is of second order for 3 < D < 4 . For 2 < D < 3 two types of the phase 
transition, first and second order, coexist. 

The dynamically generated fermion mass at low temperature and chemical potential 
disappears above certain critical temperature and chemical potential. In the following 
we shall examine the behavior of the dynamical fermion mass as a function of the tem- 
perature and chemical potential by numerically solving the gap equation. As typical 
dimensions we choose D = 2.5, 3.0, 3.5 . 

In Fig. 15 we show the behavior of the dynamical fermion mass as a function of the 
chemical potential for some fixed temperatures. We clearly observe that for D = 2.5 
the mass gap appears at the critical chemical potential as T increases while no mass gap 
shows up above D — 3 . This behavior at D = 2.5 is a reflection of the smooth transition 
from the first-order to second-order phase transition. It should be noted here that, while 
at D = 3 the phase transition is of second order, there is an exceptional case at T = 
where the mass gap exists at the critical chemical potential /i = m. In Fig. 16 we present 
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the behavior of the dynamical fermion mass as a function of the temperature for some 
fixed chemical potentials. Here for D = 2.5 the mass gap observed at low temperature 
disappears for higher temperatures while no mass gap is observed above D = 3 . 
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Fig. 15 Dynamical fermion mas 
mp^ as a function of the chemi- 
cal potential ji with temperature 
kT/m fixed at 0,0.2,0.4. 
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Fig. 16 Dynamical fermion mass 
info as a function of the tempera- 
ture T with the chemical potential 
n/m fixed at 0,0.4,0.8. 



4 CONCLUSION 

We have investigated the phase structure of four-fermion theories at finite temperature 
T and chemical potential /x in arbitrary dimensions by using the effective potential and 
gap equation in the leading order of the 1/N expansion. The theory under consideration 
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is renormalizable below 4 dimensions and give an insight into the phase structure of the 
theory in 4 dimensions. 

Starting from the theory with broken chiral symmetry at vanishing T and fi we 
calculated the renormalized effective potential for finite T and fi in the leading order of 
the 1/N expansion in arbitrary dimensions 2 < D < 4 . We found that the broken chiral 
symmetry was restored at a certain critical temperature and chemical potential. The 
phase transition from the broken phase to the symmetric phase is either of first order or 
of second order for 2 < D < 3 and is of second order for 3 < D < 4 . We found the 
boundary curve dividing the symmetric phase and broken phase in the T-\i plane. The 
dynamical mass generated in the broken phase is studied as a function of the temperature 
and chemical potential. 

At D = 2 and D = 3 formulae derived in the present investigation mostly reduce 
to the known results in the preceding works. Some results obtained in the present work 
are not known previously. On the critical curve we found analytic expressions for some 
specific points. 

Although the present work is restricted mostly to the analysis of the mathematical 
properties of the four fermion theory, we are interested in applying our results to physical 
problems in the early stage of the Universe. We will continue our work further and hope 
to publish reports on these problems. 
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APPENDIX A 



EFFECTIVE POTENTIAL FOR T = fi = 



In this appendix we present details of the calculation of the effective potential given in 
Eq. ( |2.6|) . After the Fourier transformation Eq. ( |2.5|) becomes 



— a 



(A.l) 



It should be noted that the effective potential is normalized so that Vo(0) = 0. The 
second term of the right-hand side of Eq. ( |A.1| ) is rewritten as. 



d D k , , (H 
In det ■!- 



a 



(2n) D 

l r d D k 

~2 



1 

— tr 
2 



(2ir) D i 
d D k 



In det 



— a 



+ In det 



In 



-fc 2 + a 2 
-k 2 



(A.2) 



Performing the Wick rotation k° — > we get 



1 

— tr 
2 

1 

— tr 
2 



(2^F l 1 + & 



d u K r 1 



(27T 



dx 



1 



-r i 



D 



D 



a 



(A.3) 



(2tt) d / 2 D V 2 

Inserting Eq. (|A.3|) into Eq. (|3.3|) the effective potential ( |2.6|) is obtained. 

Next we consider the two, three and four dimensional limit of the effective potential. 
Taking the two dimensional limit D — > 2, we get 



Vojcr) 



la 2 1 a 2 I , a 2 N 
^t^t + ^ j -3 + In— 



2A o- 2 47T cr s 



(A.4) 



where we use the renormalized coupling constant A defined in Eq. ( 2.10Q with 
Taking the three dimensional limit D — > 3, we get 



Vo(o) 



2\a 2 



3ao 3 , 



(A.5) 
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If we take the four dimensional limit D — > 4, we find 



2 



< 2A a 



(4< 



2 



6 ( 7 + ln2vr + - ) — 7 + ln27r + - - In— — 



e 4 3/ a 2 \ e 4 2 (Jq 2 / (J 



-0(e 4 ) , (A.6) 



where 



e 4 = . (A.7) 

To consider the theory in the case of D = 4 as a low energy effective theory stemming 
from the more fundamental theory we calculate the effective potential using the cut- 
off regularization in the case of D = 4 and compare it with the one obtained by the 
dimensional regularization. We use the renormalization condition (|2.7|) and calculate the 
effective potential ( |3.3|) by the cut-off regularization 



Vb(cr) 1 a 2 1 

(t 4 " 2A(T 2 (4tt) 















V a o 2 







A 2 1 , ff 2 \ a 1 



where A is the cut-off parameter, which is assumed to be larger than a and a®. This 
result show us that there is a correspondence between e 4 and A in the leading order of 
the 1/N expansion, 

1 A 2 

7 + ln2vr + 1 <-> In — . (A.9) 

e 4 cr 

APPENDIX B 

EFFECTIVE POTENTIAL AT FINITE TEMPERATURE AND CHEMICAL POTEN- 
TIAL 

We shall show two types of representations for the effective potential and the gap equa- 
tion at finite temperature and density in this appendix. One represented by a momen- 
tum integration is convenient for numerical calculations and the other represented by a 
summation is convenient for analytical calculations. In the Feynman functional-integral 
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formalism, thermal Green function is obtained by the following way. 



(a | e-^-^ I a) = (a, -i(3 | a, 0) 



C I [#][#] exp (ij 1 dt J d D ~ 1 x(£ + fiNU , 



where the number operator N is expressed 

N = J dP^x^i) . 

Following Eq. 
sity. 



(B.l) 



(B.2) 



we obtain the generating functional at finite temperature and den- 

(B.3) 

(B.4) 



Z^[0] = J [#][#] exp (ijf ^ dt I d D - l x(C + iiN) 
Thus the effective potential ( p.3[ ) is modified as 

V{a 



1 1 °° 

——a 2 — — 

2A P 



d D 1 k {-iw n j° - ikcf - a + //7 ) 

In det- 



(27r) D 1 (— iWrcf — ikiY + ^7°) 

where the roman index % is taken over the space components (i = 1 ~ 3). We calculate 
the integrand in the second term. 



In det- 



-i(w n + ijiyf - ikff - a] 

[-i{w n + 2/1)7° - ihrf] 

r , s 
tr / ds- — — — 

Jo [w n + + kik' 1 + s z 

1, , (w n + i/j,) 2 + kit + a 2 
-tr m — — — — — . 

2 (w n + iii) 1 + kik 1 



tr / ds 



1 



o i(w n + z/i)7° + ikff + s 



(B.5) 



To obtain the second line we use the well-known property that the gamma matrix is 
traceless. Hence we get the effective potential at finite temperature and density. 



1 1 oo 



d D ~ x k , (tu n + in) 2 + k^ + a 2 
-tr m- 



(2tt 



,£>-! 



(B.6) 



w n + in) 2 + /cjA; 1 

Next we explain how to deal with the summation and integration in the effective potential 
( B.6 ). We can perform a summation in the following way. We divide the summation into 
two parts. 

^ , \{2n+l ^ 2 
2^ ln 







+ ill 
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n=0 

oo 

n=0 



2n + 1 



2n+ 1 



n=0 



-2n - 1 







2n+l 



2 



+ ifi] + b 2 



+ (b + fif 



Using the formula 



£ In (a\2n + l) 2 + b 2 ) = In cosh— , 



n=0 



2a 



we can perform the summation and get 



_2-D/2-l_ 



1 r d u ~ Y k 



(3 J (2tt) 



D-l 



(27 

In r-^ ; h In- 



1 _|_ e -/3(|fc|+M) 



1 + e -/3(|fc|-M) 



We can also integrate over space components in Eq. ( |B.6| 
d D ~ l k (w n + i/J,) 2 + kik* + a 2 



(2tt) d - 1 (w n + iii) 2 + 



d D - x k r l 



dx 



(2ir) D - 1 Jo o 2 x + (w n + ifi) 2 + kik 1 
1 - D 



(4 7r )(I3-l)/2- 



X 



Hence we get 



1 2 1 



/?V / 2(2tt)( d - 1 )/ 2 



1 — D 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



x x: [{K+./i) 2 + ( x 2 }^- i )/ 2 -{K+^) 2 }( D - 1 )/ 2 ] . (B.n) 



From the two representations of the effective potential ( |B.9| ) and ( |B.11| ) we can write the 
gap equation in the two forms. 
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2 D/2- 



d D - l k 
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£ {K + ^) 2 + a 2 }( D - 3 )/ 2 = 0. (B.13) 
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